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Abstract 



The possibility of formation for a bound state of a t quark and a lighter 
one is investigated using potential model predictions and heavy quark ef- 
fective theory approach. Resulting estimates for the 1S-2S splitting of the 
energy levels are compared to the total top decay width T t . As for the 
case of toponium, our conclusions show that the probability of formation for 
q | T-mesons is negligibly small due to the high top mass value. 
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1 Introduction 



It is conventional wisdom JIJ that the top quark has no probability of meson 
formation of any kind because the t quark decays in a way that is too fast 
to allow for a single orbit of the bound state. First calculations based on a 
relatively light top quark || have shown that the observability of ti bound 
states would be possible for a narrow window on m t values. Subsequent 
estimates || |4[] have demonstrated that there is only a small probability 
for creation of a tt bound state for higher top masses. From these results 
it could be inferred that even the formation possibility for a meson made 
out of a single t quark and a lighter one is small. In this note we discuss 
quantitatively the last point in light of the present high values for top mass 
[|, | . For this purpose we shall make use of calculations from potential 
models in QCD in both nonrelativistic and relativistic theories. The choice 
of the potentials is driven by considerations on the proper QCD scale for this 
problem. We will also give some results taken from heavy quark effective 
theory, based on an estimate for the inertia parameter A . The latter method 
is particularly suitable for mesons containing a very light quark. 

These results will be compared to the decay probability of the top quark, 
whose value is determined by m t , in order to give quantitative answers on 
the possibility of creation for such "superheavy" mesons. 

The paper is organized as follows: we review the toponium case in Sect. 
2. The potential model approach for the T-mesons, both nonrelativistic 
and relativistic, is presented in Sect. 3. The heavy quark effective theory 
discussion for these mesons is in Sect. 4. Sect. 5 is devoted to the discussion 
and conclusions. 

2 Toponium bound states revisited 

The recent discovery of the top quark by CDF |5[] and D0 collaboration |J, 
and subsequent measurements |7| and || have given large values for the top 
mass: 

m t = 175.6 ± 9.3 GeV (CDF) ; m t = 169 ± 11 GeV (D0) (1) 

Previous top mass values of around 130 GeV already excluded the possibility 
of formation of toponium bound states || |[ . 
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This fact can be understood by comparing the time period of the would- 
be bound state with the decay probability of the top quark. For t mass above 
the Wb threshold, the top decays into a real W and a b quark || [TTJ. The 
single quark decay width of the top quark, with one-loop QCD correction 
and neglecting terms of order (mb/m t ) 2 , is given by [|ll|] : 
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The width values increase from 438 MeV (for m t = 120 GeV) to 1.02 GeV 
(m* = 160 GeV), and 2.22 GeV for m t = 200 GeV . With such a large width, 
the lifetime of T-mesons and toponium are dominated by single quark decays. 

For a bound state, the typical formation time of a hadron is characterized 
by a revolution time driven by strong interactions. Thus, no bound states 
exist, if the revolution time, tn = 2ixr/v, is larger than the lifetime of the 
rotating quarks, r ft = 1/T ft [g] . 

The value for the width of the toponium system is two times the width of 
the single top quark, since each one could decay in an independent manner: 
T t j = 2T t ■ To discuss the bound state, let us begin with the toponium case 
through a Coulombic two-body potential: 

V(r) = - 



a. 



(3) 



for which analytic solutions exist. Here should also be considered corrections 
from Higgs boson exchange Yukawa type forces [12, [L3| . For a top quark 
with mass less than w 200 GeV these corrections are small, with an attractive 
potential for the quark-antiquark singlet state, and they amount to no more 
than 10% of the Coulomb term. 

We shall use the two-loop expression for a s [[HJ , 



a s (Q 2 



47T 



2fl log [log [Q 2 /Ajts]} 
(51 log [Q 2 /A 



MSI 



with j3i 



1 A) log [Q 2 /^jts\ I 

o = ii - |n/ , Pi = 51 — Y n f evaluated at a fixed scale Q 2 
where is the Bohr radius 
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In figure ([]]) we show the behavior of the Bohr radius as a function of the 
quark mass for a quarkonia bound state under the effect of a Coulombic 
interaction. The energy is given by 



En = (6) 



8 \ia\ 

9 ~n 2 

where [i is the reduced mass. 

To estimate , we shall make use of the virial theorem, which, for the 
Coulombic potential states that (T) = —\{V) (T and V are respectively 
the kinetic and potential energies). From the energy expression (|6|), we have 
(v 2 ) = 8/9 a 2 , which leads to t R = \/2tx 9/(4a 2 s m t ). As T t j « 1crr? t (c is a 
constant) , this existence criterion allows for toponium formation if r t j > £r, 



namely for m t < y a/2/ (9tt) a s / a/c . 

Another slightly different criterion |4j] states that the formation of a 
hadron can occur only if the level splitting between the lowest lying levels of 
the bound state, which depends upon the strength of the strong force between 
the quarks and their relative distance IJ , is larger than the natural width of 
the state. In this case, we should have a bound state if AE 2 s-is > ^tt ■ 
From @ we read the splitting of the levels, AE2S-1S — 4:/3m t a 2 . In this 
case, the estimate for toponium existence is for values of top quark mass such 
that m t < 1/a/6q; s /a/c . 

A comparison with the previous method shows that the mass bound is 



larger for a factor y37t/(2a/2) , that is approximately two times larger, and 
we conclude that the energy level splitting criterion is less stringent than the 
comparison of the revolution time. We decide therefore to employ the looser 
condition: if the formation is prohibited by the AE2S-1S criterion, then it is 
certainly prohibited also from the other method. 

In figure (|2]), we show a comparison of the energy splitting and the to- 
ponium width T t j , as a function of the top mass. We shall use the value of 
Ajjs for which a s (M z ) = 0.118 and we use for a s the scale (|5|). From 



this figure, it is possible to see that for most recent top masses and the av- 
erage value for the top mass, m t = 172.9 ± 7.0 GeV, toponium bound state 
formation is rather unlikely. 
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3 T— Mesons from potential models 



For the T-mesons, the less stringent criterion of level splittings will be used. 
Notice that from f3| the comparison of the Coulombic potential results with 
some other models gave substantially similar results, thus confirming the 
validity of this approach. 

The problem for the T-meson, unlike toponium, is difficult due to the 
identification of the proper QCD scale for this problem in the Coulombic 
potential. The interquark force cannot clearly be set simply at the top mass 
scale, but it has to include both the information of the top and the light 
quark mass scale. In order to achieve it, we use for the scale the inverse of 
the Bohr radius (|5|) calculated in terms of the reduced mass of the T-meson, 
which in turn, because of the heavy top, is approximately the light quark 
mass. In addition, we should also include calculations using other potentials, 



such as Martin's [16 



V(R) = -8.064 + 6.8678 r ai (7) 
(the units are in GeV), and the model of Grant, Rosner and Rynes [TTj 



V(r) = -{r a -l) + c (8) 

a 

These have the property of being independent of QCD coupling a s , an 
important consideration since the energy scale considered is low, of the order 
of few fractions of GeV. For quark masses we use the constituent one, whose 
approximate values are m;, = 5.0, m c = 1.5, m s = 0.5, m u = = 0.3 (in 
GeV). In table (1), we present a table of energy splitting values, for the bound 
states from the potential models. For sake of comparison we also include the 
prediction of the Coulombic model and a s values. As for toponium, we shall 



use Ajmy such that a s (Mz) = 0.118 GeV [[0|, and the scale from 
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Energy splitting values AE for potential models (GeV) 



quark 


Martin 


Rosner 


Coulomb 


a s 


b 


0.56 


0.65 


0.29 


0.30 


c 


0.60 


0.61 


0.20 


0.45 


s 


0.31 


0.41 


0.20 


0.78 


u, d 


0.08 


0.02 


0.24 


1.09 



Table 1: values of AE from different potential models (PM). We show 
AE as a function of the light quark mass for the potentials given in for- 
mulae^), (g) and (@). 



We see that AE increases with the mass of the light quark, except for 
the Coulombic case, whose high a s values for light quarks make perturbative 
calculations unreliable. 

In the following, we investigate the bound state system by means of a 
model that includes relativistic corrections. The t quark inside the tq system 
moves nonrelativistically ( v 2 /c 2 ~ 0.01 for m t « 173 GeV ). In order to 
show whether the lighter quarks need a relativistic treatment we employ the 
Salpeter equation: 



-V 2 + m 2 t 



-V 2 + m 2 + V(r) 



4> = Eip 



(9) 



with potentials from the nonrelativistic models. Instead of dealing with 
the mathematical difficulties of finding the eigenvalues of flj^) because of the 
square root operator, we shall apply the Rayleigh-Ritz variational method 
with suitable trial functions. This procedure has already been successfully 
applied for the B and D meson decay constants |TB[ . By means of a Fourier 
transform on fl9|), we have 



p 2 + m 2 + Jp 2 + m 2 + V(r) 



ip = Ei/j 



(10) 



Writing explicitly the dependence of ip upon an extremum parameter £ , 
we find that the energy of the state is given by minimizing the expectation 
value of H in (Bl) 
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(H) = mO\HmO) = E(0; ^ = Ofor£ = £ (11) 

The energy levels depend on the number of nodes of the trial function: no 
nodes for the ground state, and one node for the first excited level. We have 
performed the calculation for two kinds of trial wavef unctions: hydrogen- 
type wavefunction, coming from a Coulombic potential, and Gaussian wave- 
function, from the harmonic oscillator, thus mimicking the short and the long 
range behavior of the interquark force respectively. 
For the IS* hydrogen-like function, we choose 



and its Fourier transform, which reads: 



2 \/2 a 3 / 2 

^ = -IT n f 2 2,2 ( 13 ) 

7T (1 + a 2 p 2 ) 



a being the variational parameter. The 2S function is: 

J_ a -3/2 ( x _ £. 
V8^ V 2a 



^( r ) = ^ a " 3 / 2 (l-— )e- r ' 2a (11) 



and 



fe) = 15 a ^ M^_I> (15) 

7T (l + a^p 2 ) 



respectively. For the Gaussian wavefunction, the IS* is: 



^(r) = UA e-^ 2 ' 2 (16) 



and 



(v^) 3/2 

The 2S function is given by: 



$(p) = — — ™ e^ 2 / 2 - 2 (17) 
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and 



(5v^F-8)7r 



3/2 



ru) e 



V^/2 



(18) 



$(p) 




(50F-8)tt 



-3/2 -pW 



(19) 



where \i is a different minimization parameter from the previous one. L a n {x) is 
the generalized Laguerre polynomial, obeying the differential equation xy" + 
(a+1— x)y'+ny = and the orthogonality relation / °° L a m {x}L a Jyx)x a e~ x dx = 
for m 7^ n (for a = , we retrieve the usual Laguerre polynomials). 
Technically, one splits H = T + V and then computes the average 



E(0 = $(p)\T(p)$(p) + (V(r)|V(r)Mr)) 



4tt !™dpp 2 ^(p) 



p 2 + mf + Jp 2 + m- 



i)(p) + 47r / °° cir r 2 ^(r^r)^?") 

(20) 



Below, we present a comparative table with results from both the potential 
models and two kinds of wavefunctions. 



AE for potential and relativistic models (GeV) 





PM 


RPM, Coulomb 


RPM, Gauss 


quark 


Martin 


Rosner 


Martin 


Rosner 


Martin 


Rosner 


b 


0.56 


0.65 


0.47 


0.57 


0.53 


0.61 


c 


0.60 


0.61 


0.45 


0.39 


0.55 


0.35 


s 


0.31 


0.41 


0.41 


0.09 


0.44 


0.00 


u, d 


0.08 


0.02 


0.38 


0.01 


0.34 


0.00 



Table 2: Comparison of AE results from potential models (PM) and 
relativistic potential models (RPM) for potentials (0) and dj|) . "Coulomb" 
and "Gauss" labels indicate the two different trial wavefunctions. 
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As we should expect, major differences between the nonrelativistic and 
relativistic potential models arise for the lighter (and hence faster) quarks. 
Also, the Martin potential is less sensitive to mass changes. We also no- 
tice that there are no significant differences between different kinds of trial 
wavef unctions; confirming thus the reliability of this approach. 



4 Heavy Quark Theory approach 

For performing these calculations, we shall now use a different method bor- 
rowed from heavy quark effective theory (HQET) . This has the advantage 
of giving results at a scale (u and d quarks) for which some potential mod- 
els predictions may not be reliable enough. The Coulombic model is one of 
these, because a s is of order 1 in this case. 

In this model, a meson (hadron) containing a single heavy quark (mq ^> 
Ajjs) is considered. The heavy quark's momentum can be written as 

p Q = m Q -v + k (21) 

k is the "residual" momentum, which measures the degree to which the quark 
is off-shell [I9|; v is the velocity satisfying v 2 = 1 . The quark Q exchanges 



only small momenta with the rest of the hadron, so it is essentially on shell, 
Pq = rriQ . Q behaves like a static electric and chromomagnetic field source. 
The properties of the light degrees of freedom do not depend upon the flavor 
and mass of the heavy field. Since the top quark is very heavy, it plays the 
role of the heavy quark and we could apply HQET to give some information 
on the energy of the bound state. 

The mass M of the meson is here expressed as 



M = m Q + A + o(^^J (22) 

where A is a positive contribution to M . The "inertia" parameter A has no 
dependence on the heavy degrees of freedom: 

A = A„ = lim (M - m ) (23) 

mq-»oo 

On the other hand, from the potential models (PM) it is possible to write, 
for the meson mass, 
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M = m Q +m q -E b (24) 

where m q is the mass of the light quark, E b (> 0) is the binding energy [[| . 
Comparing fl22|) and (p4|), we arrive at the identification 



A g = m g - E b (25) 

neglecting terms of order 0(l/mg) , an operation valid for m t ~ 173 GeV . 
We find therefore that the quantity A q gives substantially the energy of the 
light degrees of freedom^. In order to compute AE, we need some estimates 
for A„ . 



A rigorous lower bound on the value of A g has been derived |22j, and 
leads to the values (in GeV) 

A~ q (Qu) > 0.057 , A~ q (Qd) > 0.076 , ~A~ q (Qs) > 0.343 (26) 

for a meson formed by a heavy quark Q and an u, d and s quark respectively. 
This lower bound has been obtained using the Euclidean path integral for- 
mulation of QCD. The Cauchy-Schwarz inequality has been used to derive 
inequalities among Euclidean correlation functions, obtaining m(QTq)—m q > 
1/2 m'(qi , -f5q) . The first term of the LHS is the heavy meson mass, while 
the RHS is the mass of the meson created by the light degrees of freedom, 
neglecting annihilation diagrams. 

Equation (|^) gives us the binding energy of the lowest state, while we 
need the difference in energy. From the above 

m q -E b = A q > lq (27) 

where 7 g is the lower bound on A q . Moreover, we have E b > AE 2 s-is > an d 
we obtain from the lower bound on A q , an upper bound for the quantity 
AE 2 s-is '■ 

AE 2S _ 1S <m q - lq (28) 



From equation (P4|) and constituent mass values for the light quarks, we 



obtain the following results: 

AE 2S -is(tu) < 0.243 , AE 2S -is(td) < 0.224 , AE 2S -i S (ts) < 0.257 (29) 
2 A similar equivalence has been obtained in pit] . 
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with m u = ma = 0.3 GeV ; m s = 0.5 GeV 



The values in equation (29) have been derived using chiral perturbation 
theory, and such a method is not applicable for heavier quarks. There is still 
the bound A g > ^m'(qi^y 5 q) p2j , and therefore A g > . With this less 
stringent bound and using the aforementioned values, we have (in GeV): 

AE 2S „ ls (tc) < 1.5 , AE 2S . ls (tb) < 5.0 (30) 

assuming m c = 1.5 GeV, nib = 5.0 GeV . We notice that the estimate of the 
HQET for the AE 2 s-is for the heavier quarks is rather large due to a poor 
estimate on A g , yet it is consistent with the potential models. 



5 Discussion and conclusions 

In table (3), we present a summary of the results which have to be compared 
to the width of the T-meson. If Ai^s-is > Ttg , then there is an oppor- 
tunity of formation for the bound state. Since the quarks could decay in 
an independent manner from each other, one has T t g = T t + T q . We know 
however that T q w m 5 q for a light quark, while T t w mf , leading to T t ^> T g , 
thus T tq is essentially given by the top width. 



Summary of AE results (GeV) 



quark 


PM 


RPM 


HQET 


b 


0.29-0.65 


0.47-0.61 


0.00-5.00 


c 


0.20-0.61 


0.35 - 0.45 


0.00-1.50 


s 


0.20 - 0.41 


0.00 - 0.44 


0.00-0.26 


u, d 


0.02 - 0.24 


0.00-0.38 


0.00 - 0.24 



Table 3: Summary results of AE ranges from potential models (PM), rel- 
ativistic models (RPM) and heavy quark effective theory estimates (HQET). 



The tables (1), (2) and (3) have to be compared with the width values 
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of the top. Considering the average value m t = 172.9 ± 7.0 GeV for the top 
mass, we have from (0) (in GeV): 

T t = 1.17 (mt = 165.9), r t = 1.35 (m t = 172.9), r t = 1.55 (m t = 179.9) 

(31) 

Figure (|3|) shows the window in energy ranges with respect to the decay 
width of the T-meson, which is drawn in full line. The dashed horizontal 
lines represent the indetermination due to the error on the top mass. 

We see that the values for Tt are larger than the one predicted from both 
the PM and the RPM. It must be noticed that some of the values from 
the HQET, the one for the heavier quarks (see table (3)), seem to allow for 
T-meson formations, but it has to be stressed that the estimates for the c and 
b quarks are rather poor, since at the present time a more precise estimate 
on A Cj b is lacking. 

From a combination of these calculational techniques we may conclude 
that there is little evidence of possibility of formation for bound states of t 
quark and lighter quarks. 
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Figure 1: Bohr radius versus quark mass 
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Figure 2: Toponium width compared to energy separation of the Coulombic 
model. The stripe represents the indetermination of AE due to the value of 
a s . The vertical bars show the indetermination of the top mass. 
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Figure 3: Energy splittings from non relativistic and relativistic potential 
models for different light quark masses, are compared with the T-meson width 
(full and dashed horizontal lines). 
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